Motifs haven’t been updated in this yet ref = “total_ratio” changed to “self” when using show_burden_for_selected_features Use a newly calculated set of motif change scores to define motif SNVs, but only focusing on motifs that more relevant Here motifs are defined as relevant if they are 1) in top 50% brain expressed genes, and 2) are enriched in at least one gain-enriched brain module Partition enhancers based on distance to TSS based on Yanyu’s pipeline. Basically, her pipeline takes intersections between enhancers and distance ranges, and remove bases of enhancers that are not overlap.

non ASD genes are the bottom 20% not plotting mutation is 5’UTR or 3’UTR added one row for showing the union of enhancers and promoters Only have Fantom, Noonan and Whole genome Look at more genesets.

getwd()
## [1] "/media/yuwen/F/ASD/analysis"
knitr::opts_chunk$set(warning=FALSE, message=FALSE)

read in coding tada data

## [1] "/media/yuwen/F/ASD/analysis"

Whole genome

Noonan enhancers

Fantom fetal brain

Brain roadmap adult brain

Brain roadmap + Noonan brain adult brain

plot_summary = plot_summary[!is.na(plot_summary$gene_list),]
plot_summary = plot_summary[!duplicated(plot_summary),]
#scaling factor for plotting, which is the burden of baseline. Used for 160229 data only
scaling = 0.73